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We find a one-parameter family of variables which recast the 3+1 Einstein equations into first- 
order symmetric-hyperbolic form for any fixed choice of gauge. Hyperbolicity considerations lead 
us to a redefinition of the lapse in terms of an arbitrary factor times a power of the determinant of 
the 3-metric; under certain assumptions, the exponent can be chosen arbitrarily, but positive, with 
no implication of gauge-fixing. 
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The issue of setting up a well-posed initial-value formulation for general relativity has been studied with the help 
of varied strategies, including special gauges and higher-order formulations 0. Recently, a renewed interest Q in 
the problem has arisen, in connection with the numerical evolution of the Einstein equations away from an initial 
hypersurface. Although the relevance of a manifestly hyperbolic formulation to the numerical integration of the 
Einstein equations is not yet clear, it is believed that a code tailored in a hyperbolic formulation would share properties 
of the exact system; namely, it would guarantee uniqueness and stability of solutions evolved from proper initial data. 
However, technical issues, associated with the discretization of the equations and the precision of the approximation, 
which may concern numerical stability, are not necessarily ruled out by a pure theoretical hyperbolic development. 

In regards to the manifest hyperbolicity of the Einstein equations, the relevance of gauge choices has long been 
a question open to consideration. The gauge freedom of 3+1 general relativity is embodied by the lapse function 
and shift vector, which are completely arbitrary since their evolution is not determined by the theory. In general, a 
theory expressed in terms of equations on fields which admit gauge freedom may not admit a well-posed formulation 
unless improper gauge choices are ruled out, or the true gauge-invariant variables of the theory are found. A typical 
example is Maxwell's theory; on the one hand, it admits a hyperbolic formulation at fixed gauge and in terms of 
gauge invariant variables as well; on the other hand, anomalous gauges can be found for which the resulting system 
does not have a well-posed initial-value formulation. 

Our intention is to give an explicit argument to rewrite the 3+1 Einstein equations into a manifestly well-posed 
form, without the need of resorting to a choice of gauge. It has certainly been known that general relativity in special 
gauges can be set in symmetric hyperbolic form ||. Furthermore, it has recently been shown that, for certain special 
first-order variables, general relativity admits a symmetric hyperbolic formulation for arbitrary fixed gauge QJ^]. Here 
we extend the existing results by showing that, under certain assumptions, there is a one-parameter family of new 
first-order variables for general relativity which satisfy first-order symmetric hyperbolic evolution for arbitrary but 
fixed choice of gauge. 

Several concepts of hyperbolicity (e.g., strict, strong) can assert the well-posedness of a system of PDE's. Among 
all different concepts, symmetric hyperbolicity is especially appealing, for the reason that most interesting physical 
systems admit a formulation of this type |(|. Symmetric hyperbolicity is based on the symmetry properties of the 
differential operator 0; therefore, multiple eigenvalues, which usually occur due to the presence of symmetries, play 
no role, as opposed to the case of other types of hyperbolicity. The reason for the well-posedness in the symmetric- 
hyperbolic case is that the energy norm (an integral expression in terms of the fields) at later times can still be seen 
to be bounded by the norm at the initial time, because of cancellation of terms under integration by parts. The 
symmetry of the differential operator in the evolution equations guarantees the cancellation. 
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In the following, we set up the problem of general relativity in the 3+1 formulation due to Arnowitt, Deser and 
Misner in a non-canonical (though widely used) choice of variables; i.e., the intrinsic metric and the extrinsic curvature 
of the spatial hypersurfaces. We then redefine the variables in order to reduce the system to first order; the redefined 
variables depend on a set of parameters to be fixed by hyperbolicity considerations. Finally, we use the argument of 
the cancellation of terms under integration by parts in the energy norm to determine the parameters. In the process, 
we find that the lapse function must be redefined in terms of the determinant of the 3-metric, without any loss of 
gauge freedom. 

In order to fix notation, we summarize some necessary points of the 3+1 formulation. The 3+1 splitting of the fully 
4-dimensional formalism consists of a spacelike foliation by the level surfaces of a function t{x a ). The unit normal 
form is n a = — iVV i, where N is the lapse function. The unit normal vector is given by n a = -jj{t a - N a ), where 
7V a is the shift vector. The metric g ab induces a 3-metric on the spatial surfaces, by h ab — g ab + n a n b . In a coordinate 
system {x°,x 1 }, i = 1,2,3, adapted to the surfaces (such that t — x°), the induced metric h ab reduces to h % K The 
extrinsic curvature of the 3-surfaces is defined by K ab = ^£ n h ab , and is also a 3-tensor. 

The equations for the evolution of the intrinsic contravariant metric h %3 and the extrinsic curvature K lJ can be 
taken as (from |s) , with the notation of Chapter 10 and Appendix E of Q) 

h ij = 2NK ij - D i N j -D 3 N l (1) 
k" = A ( //'•' - 2K lk K 3 k + K ij K 

- -h ij (S - /?))) - D l D 3 N 
+N k D k K lj - K lk D k N j - K ]k D k N l (2) 

The notation ( ' ) stands for £f or simply d/dt. Indices k, . . . are raised and lowered with the 3-metric h lJ ] the 
operator Di is the covariant derivative with respect to the 3-metric hij |jj . For any 3-tensor L™ , the notation U stands 
for its trace with respect to the 3-metric: U = U k k . The matter tensor S lJ is the projection of the 4-dimensional 
stress-energy tensor T ab into the spatial hypersurfaces, and p is the projection of T ab in the direction normal to the 
surfaces. The results derived in this Letter do not depend strongly on the particular matter source, but hold for any 
sources that admit a first-order symmetric-hyperbolic formulation on their own. 
Equations ([j]) and for the fields , K % i ) are supplemented by the constraints 

C := ^(R + K 2 ~ K zj K i: >) - k P ^0, (3) 
C := DjK^ - D l K - nf = 0, (4) 

where J 1 is the mixed projection of T ab onto the hypersurface and the normal. If C and C 1 can be shown to be 
conserved as a consequence of (Q) and (||), then the constraints only need to be imposed on an initial hypersurface. 
This will be our point of view in the following. 

Introduction of the parameters. In order to show that there exists a one-parameter family of variables for which 
general relativity takes a first-order symmetric-hyperbolic form, we first introduce a set of four parameters, a, (3, 7 and 
e; we eventually require the parameters to satisfy a set of three algebraic conditions that guarantee the hyperbolicity. 

Two parameters, a and (3, are used to redefine variables as follows 

APh = l(h^, k + ah^h rs h rs , k ), (5) 
pij = Rij + p h ij K. (6) 

The definition (||) reduces the Einstein equations to first order. Note that the variable M % \ represents the spatial 
derivative of the densitized 3-metric, k = ^h a (h~ a h %: i}, k , where h is the determinant of hij. Equations (||) and 
(^) can be inverted into 

h», k =2(M*S k - -^-h ij M k ), (7) 
da + 1 

K i3 = pi3 P—h l3 P, (8) 

3/3+1 ' v ; 

with the notation M k = h^M n k . 
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A third parameter, 7, is introduced in the evolution equations to allow for a combination of ([!]) and (|J) with the 
constraints (|) and (Ej). In this way, the principal part of the evolution equations ([j]) and (||) can be modified. The 
constraints (|3) and (|j) will be assumed to be conserved by the resulting equations. Since 7 plays a crucial role in the 
hyperbolicity of the system, in the following wc point out its exact place in the evolution equations. 

The evolution equation for M lj k can be obtained by, first, taking a space derivative d/dx k of Eq. ([!]), and then 
tracing and combining the resultant equation according to the definition (||). We also add the vector constraint C l 
with an appropriate -uniquely determined- factor, obtaining the following: 

M t3 u = l((h ij ),k +a(h ij h rs + h^h rs )h rs , k 

+aWh rs (h rs ) lk ) - N5%Cfi . (9) 

If Eq.(Q) is used in the right-hand side of (^|) to eliminate time derivatives in favor of space derivatives of the new fields, 
the right-hand side becomes a combination of the fields , M y t, P y , lapse, shift and sources; first-space derivatives 
of M lJ fc,P tJ , lapse and shift; and second-space derivatives of the shift. This equation is shown explicitly below [Eq. 
(|l7|)], correct to principal terms. 

The evolution equation for P y is obtained directly from Eq.fQ), by the appropriate combination with its trace, as 
prescribed by the definition (H). We also add the scalar constraint C with a suitable factor: 



pn = %v + p{h %3 K + h l: >h rs K rs + h i: >h rs K r 

+2N-yh ij C. (10) 



When Eqs. (|lj) and (|2|) are substituted appropriately in (|10|), the right-hand side becomes a combination of the 
fields /i'^M y t ) P' J ', lapse, shift and sources; first-space derivatives of M y 'fc,P* J ', lapse and shift; and second-space 
derivatives of lapse. This equation is shown explicitly below [Eq. (fL6|)], correct to principal terms. 
The fourth parameter, e, is introduced in order to redefine the lapse N by 

N = h- i3a+1 ^/ 2 Q , (11) 

for an arbitrary function Q. The lapse is thus redefined without loss of generality; the gauge freedom is transferred 
to Q, and the parameter e remains to be specified. Notice that 

^ = eA4 + (lnQ), fe . (12) 

Since second derivatives of the lapse appear in (|l^), this redefinition allows for a modification of the principal terms 



in (10) 



Hyperbolicity imposed on the system. Wc define the energy norm of the system at time t as 

E(t) = //'•'//,., + /'"/',., + M^ k M l} k (13) 

where the integration is performed on the surface S defined by t = const. The spatial symmetry of the system is 
guaranteed if the principal terms in the time derivative of the energy (Tc| l can be combined into total divergences, 
since in this case their contribution to the energy estimates would vanish. 
The time-derivative of the energy, correct to principal terms, is 

E(t) = f Wtiij + /"••/',., + M t] k M l3 k . (14) 



The evolution equations (|l|), (|9j) and (10) can be used to trade time derivatives for space derivatives in (|l4j). If the 
principal terms can be eliminated under integration by parts, then the system becomes hyperbolic. 

In the following we write the principal terms of the evolution equations and find the conditions that are necessary 
to symmetrize the system. 



The Ricci tensor fly is needed in terms of the new fields. Recall 

p j^k j^fc I j^k pi j^k T1Z 

— 1 ij.k ~ 1 ki.j "+" 1 ij l kl ~ 1 lj L ki I 
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with r& = — ^huh k \j — ^hjih kl 7i + ^h kl h ir hj s h rs ,i . In terms of M tJ 'k, the connection r*. takes the form 

2a 



2h l{i M kl }) + h kl h ir h 3S M rs i -r 3^i«(i-«i) 



3a+l 



(15) 



The principal part of ( [To| ) is then 



py = N k P l \ k + N(h kl M ij k ,i - 2h l{l M j)k hk 

3a + l 3a+l 
-2/3h ij M kl k ,^ - h lk W l N, kl ~f3h^h kL N, kl 

+ 2N 1 W ( - M H M + £±±-h kl M k , t ) . 
\ 3a + l / 



(16) 



The principal terms of (0) are the following: 



M %3 1 



N l M^ kll + N(P«, k + 



a — f3 



WP lk 



2N5 [ k l P :i)l ,i + 2N 



3/3+1 



P+ 1 5 (iu.vi 
3/3+1 k 



(17) 



In view of (^6|) and (^7|), the cancellation under integration by parts in ( |l4] ) takes place if the following algebraic 
conditions are imposed on the parameters a, /3, 7 and e: 



2a + 1 
3a + 1 
+ 1 



-e = 0, 



3/3+1 

2(/3 + 7 )(a + l)-a a - /3 



3a + 1 



3/3 + 1 



+ /3 + 7 = 0, 
/3e = 0. 



(18a) 
(18b) 
(18c) 



Condition ( 18a ) has the effect of the cancellation of the fourth and seventh terms in (16), even before their contribution 
to the energy is considered. This is done in this way, because the fourth term in (|l6) has no symmetric counterpart in 
( i7|) with respect to its contribution to the energy, and, therefore, needs to be eliminated from the system. Condition 
(18b ) guarantees the cancellation, under integration by parts, of the sixth and ninth terms in (16), together with 
their symmetric counterpart, i.e. the fifth term in (|l7|). Lastly, condition (18c) guarantees the symmetry of the fifth, 
eighth and tenth terms in ( |l6| ) with the third term in (|l7j), which subsequently make no contribution to E. 

With the assumptions that Q > 0, that h 1 ^ is positive definite, that the algebraic conditions ([ish are met by the 
four parameters a,/3,7 and e, and that the constraints C and C % are conserved, the fields M y fe,P y ) satisfy a 
symmetric hyperbolic system of PDE's, namely Eqs. (]l|), (||) and (10), with the initial data constrained by (^), (^) 
and (§). 

Notice that the conditions ( |l8| ) leave free one of the four parameters. Any one of the parameters can be chosen 
freely, within a real range that allows for real values for the remaining three parameters as solutions of fll8|). For 
instance, if a is considered as the free parameter, then a can take values in ( — 00, —1/2), while /3 must be chosen as 
a root of the following quadratic equation: 



3/3 2 + 2/3 - 



(3a + l)(a + l) + l 
(2a + 1) 



= 



(19) 



A most interesting choice of variables is P iJ = K 1 ^ — h l ^K (proportional to the canonical ADM momentum), or 
/3 = —1. This choice of /3, in turn, fixes 7 = 1, a = — 1 and e = 1/2 . The redefinition of lapse becomes N = QVh, and 
the number of terms in the principal parts in Eqs. ( |l6| ) and (17) reduces considerably. Regarding the propagation 
of the constraints, for 7 = 1 it can be shown that the Bianchi equations imply a homogeneous symmetric hyperbolic 
evolution system for C and C 1 . It follows that the constraints are conserved. This case was explored earlier by the 
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authors, and has been found suitable for the development of a smooth newtonian limit Q if certain gauge choices are 
imposed in addition to the well-posed formulation. Most remarkably, for 7^1 the evolution of the constraints is not 
symmetric hyperbolic nor strictly hyperbolic, and the validity of the assumption of the conservation of the constraints 
must be studied carefully. The details will soon appear elsewhere. 

Eq. (18a) shows that the exponent of \fh in the redefinition of the lapse ( pd| ) is equal to —(2a + 1), being thus 
any positive real number, but never zero. Thus, it is not possible to have a set of variables of the form (|^)-(|^) with 
symmetric hyperbolic evolution without relating the true lapse N to the 3-metric. 

The system ([[]), (g) and (0) has a non trivial set of characteristics Using the notation £ Q = (v,£i) where £i 
has unit norm with respect to h v , it is immediate to see that £ a is characteristic if i a £ Q = (for any a). Furthermore, 
by essentially the same arguments as in J5|, it can be shown that covectors £ a satisfying either £ a £ a = or n a £ a = 
are also characteristic (for any a). There are no other characteristics if a takes the value —1, as in Q. Therefore, if 
a = — 1, the characteristics are null (with speeds N 1 ^ ± N), or timelike and either tangent to t a (with zero speed) or 
tangent to the normal direction n a (with speed N 1 ^). However, if a 7^ — 1 the system may have other characteristics, 
in addition to these, with speeds that may depend (non-trivially) on the choice of a. The details will also appear 
elsewhere. 

In this work the choice of N l and Q is arbitrary but given. The gauge must be specified in order to integrate the 
equations. On the other hand, the hyperbolicity of the system holds independently of the choice of gauge. The fact 
that gauge-fixing is required is not troublesome; the variables themselves are not gauge-invariant fields. 

In order to avoid confusion, we point out that by fixing a gauge we understand a non-dynamical specification of 
N l and Q as a priori known functions for all time, independent of the evolution of the new fields. In this way, N l 
and Q act as known sources. If the gauge were specified dynamically (explicitly, or implicitly via an equation) as a 
function of the new fields, then the principal part of the evolution equations would be modified. The results proven 
here do not guarantee the well-posed evolution of such a choice of gauge; in fact, it is not hard to see that, in general, 
the well-posedness would be hampered. 

The system shown here shares, with other hyperbolic formulations, the property of manifest first-order flux- 
conservative form, which makes it suitable for the application of general numerical integration techniques ]l2| (no 
numerical applications of this formalism have been investigated as of now). Aside from that, we find it very appealing 
for its remarkable simplicity and clarity. 
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